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Q Abstract 

A new dynamic state characterized by (2m; + l)n static phase kink with integers {mi} is proposed recently in a stack of inductively 
f^coupled Josephson junctions. In the present paper, the stability of the phase kink state is investigated against many perturbations 
^ and it is shown that the kink state is stable. It is also discussed that the suppression of the amplitude of superconducting order 



< 



parameter caused by the kink is weak. 
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1. Introduction 

The experimental breakthrough fll |2l on the terahertz radi- 
'ation from a mesa of Bi2Sr2CaCu208+*(BSCCO) single crys- 
tal has triggered tremendous effort to reveal the mechanics of 
emission. Soon after the experiments, it is proposed that a 
new dynamic state in a stack of inductively coupled Joseph- 
son junctions is probably responsible for the experimental 
observationsisl |4j]. The new dynamic state is characterized by 
(2m; H- 1 )7r phase kinks that are stacked along the c direction and 
are localized at the nodes of electric field, where {m/) are inte- 
■gers. To date, this kink state is the only solution that can suc- 
cessfully explain the key observations in the experiments 
first, the radiation frequency and voltage obey the ac Joseph- 
son relation; secondly, strong and coherent emission appears at 
the cavity resonances. In the present paper, we investigate the 
stability of the kink state. 

2. Results and discussion 

■ The dynamics of the gauge invariant phase difference in 
BSCCO is properly described by the inductively coupled sine- 
Gordon equations [5, 6] 



dlP, = (1 - ^A(2))[sinP; + d^P, +/3d,P, - 7ext], 



(1) 



where P; is the gauge invariant phase difference in the Zth junc- 
tion, /ext the external current, /3 the normalized conductance 
and { ~ 10^ the inductive coupling. Dimensionless units are 
adopted and their definitions are available in Ref. [7J. Due to 
the presence of significant impedance mismatch as the case in 
the experiments|llj,|2(], the boundary condition may be approxi- 
mated as d ^Pi = in the absence of an external magnetic field. 
The kink state as one solution to Eq. ([TJ can be written asllH 



where the first is the rotating phase with voltage w, the second 
term the static phase kink and the last term the plasma oscilla- 
tion with ki = 71 1 L » 1 with L being the length of junction. 
Here we have considered the voltage near the first cavity mode. 
Because of the huge inductive coupling, P^ runs sharply from 
to (2m;-i- 1 )n in the narrowed region of width 1 / ^j^\A\. Approxi- 
mating P^' as a step function, we obtain A - 4-/7r/{a/-k^-ifioj). 
In the following, we consider the region where A < 1 . 

To check the stability of the kink state Eq. (|2|, we add a 
small perturbation to the solution, f J = f / H- 0i with 6i 1. The 
kink state is stable if the perturbation dies out and is unstable 
if it increases when f — > +oo. Substituting PJ into Eq. ([U, we 
obtain the equation governing the evolution of the perturbation 



520, = (1-^A(2)){[ 
-f cos(fcix) exp(-iP])]9, + pdfii + d^Oi). 



(3) 



Piix, f) - (jjt -¥ P^(x) -H Re[-!A cos(^i;i:) exp(;W)], 



(2) 



Here we investigate the region that w » 1, A « 1 and the 
variation of 6*; is much faster than that of PJ along the c direction 
in the following calculation. The solution of Oi can be expressed 
asH 

0i(x, t) = 

Z + exp(/<yf) H- 9q- exp(-/wf)j expiiql) exp(-/Qf)- ^'^^ 

Here the complex eigenfrequency Q is assumed to be small 
Q <K 1 . The higher frequency harmonics Q + mcu) with m > 1 
in Eq. (|4]i are small in the region of w » 1 and are neglected. 
Periodic boundary condition is imposed along the c direction, 
and hence q = 2mn/N with m being an integer and the to- 
tal number of junctions. The kink state is stable if and only if 
Imfl < for all q. 

With the condition that P-^ has slow variation along the c axis, 
we have 

(1 - fAP))[exp(/P;)0,] « exp(/P;) ^ a,e, exp(iql), (5) 
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where = I +2^(1 -cosq). Furthermore, because PJ is almost 
a step function, exp(/fp « cosP' = l-20(x-L/2)where0(^) 
is the Heaviside step function. Substituting Eq. (|4]l into Eq. Q 
and using Eq. (|5]l, we have for the frequency components w + Q 
and Q. 



cos 



-(0, 



9+ 



COS f" 



COS P" — ~ 
6a — LO.O, 



±^q± 



AcO?.(klX)0a-Q?0a 



(7) 



where a small term proportional to A in the coefficient of 6q± in 
Eq. (|7]i is neglected and 



n2 = Q2 + ySiQ; ul = (+w - or -PK+U) - Q). 



(8) 



By inspecting Eq. the length scale over which 6q varies 
is of order of 1 /(O y[a^ » 1 . To solve Eq. O, we resort to 
multimodes expansion 



cos P' = ^ a„ cos(^„x), 6q± - ^ b„± cos(A:„x), 



(9) 



n=0 



«=() 



with fl„ = sin(Y)- Substituting Eq. (|9ll into Eq. ©and 
neglecting the coordinate dependence of 9q, we obtain 



OqUl - kl 



(10) 



Inserting Eq.® into Eq. Q, we obtain the equation for 6q 

d% + aqll% 

Z (^,+ + b„-) cos(k„x) - cos(kix)e'q 



parameterlH^ . It is recognized that Eq. ([T]i assumes the am- 
plitude of superconductivity to be a constant over space and 
time. To study fully the suppression of the superconductivity by 
phase kink, one needs to start from equations self-consistently 
describing the evolution of phase as well as amplitude of super- 
conductivity. However, we assume the amplitude of supercon- 
(6) ductivity is fixed, which corresponds to infinite condensation 
energy limit, and then consider the effect of phase kink. 

The non-uniform alignment of kink along the c axis induces 
in-plane magnetic field and supercurrent. The magnetic field 
is -sJlAlJ^ <s; 1 and the in-plane supercurrent is of order of 
lOOJc with Jc the critical current along the c axis. Compared 
with the critical magnetic field and critical current density. Both 
of them are too small to suppress the superconductivity sensi- 
bly. From the energetic point of view, the energy costs caused 
by phase kink is of the order of Josephson coupling, which is 
much smaller than the condensation energy in BSCCO. Thus, it 
cannot suppress appreciably the amplitude of superconducting 
order parameter which is proportional to the condensation en- 
ergy. All these facts indicate the assumption of constant super- 
conductivity amplitude is self-consistent and the effect of phase 
kink on the superconductivity can be safely neglected. This 
estimate is consistent with the established treatment on soliton 
phenomena in Josephson iunctions jl ill : although the suppres- 
sion of superconductivity amplitude has never been included in 
junction physics, calculated IV curves for soliton states agree 
well with experimental results. If one takes into account the 
suppression of superconductivity amplitude due to the phase 
kink, the total energy is reduced from the one estimated by fix- 
ing the amplitude, which corresponds to the infinite limit of 
superconductivity condensation. 
^^^^ In short, the kink state is stable and has neghgibly small ef- 
fect on the superconductivity. 
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the spectrum of small perturbations to the kink state. Integrat- 
ing Eq.lfTTTi from to L and approximating 9q in the integrand 
as a constant when it is multiplied by other spatial dependent 
functions, we obtain the dispersion relation 



«=0 



^2 - {k'„f 



[co2-(k'„ff+pW- 



Aai 



(12) 



with k'^ = k„l -yS^ <K 1. It is straightforward to check that 
ImO < when u < ki, but solution with ImO > appears 
when 0) > kl. Thus we conclude that the kink state is stable in 
the region ofu<k\. 

To study the full spectrum of stability, we resort to numer- 
ical simulation. It is found that the state with kink is stable 
against different types of distortions, such as small perturba- 
tions both in the lateral direction and stack direction, thermal 
fluctuations, weak external magnetic field[7] and modulation 
of critical currentfl. The radiation which is neglected in the 
analytical treatment is found to reduce the stability of kink state 
to a certain degree tsj]. 

Finally, we discuss the impact of the sharp phase change 
at the kink on the amplitude of superconducting order 
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